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Abstract: Rainfall-runoff relationship in arid regions is unique and challenging to study. Studies for bridging 
the hydro-meteorological knowledge gap for planning, designing and managing water resources is therefore 
vitally important. The objective of this study is to develop a method for estimating unit hydrograph at 
reasonably finer time resolutions (10-min and 1-h) which can be easily adaptable by practitioners at sub- 
catchment levels, especially when the focus area is ungauged. Observed wadi-flow at 5-min interval and 
tipping bucket rainfall measurements at 1-min interval were obtained to cover 10 major watersheds in 
Oman. The deconvolution method was applied to derive the unit hydrographs (UHs) from wadi-flow and 
excess rainfall. Key catchment characteristics such as the watershed area, length of the main wadi and the 
length to the centroid of the catchment area were derived from digital elevation model (DEM) data. The 
whole study area was then divided into 515 sub-catchments with various shapes and sizes. A strong 
relationship was found between the wadi length and the length to the centroid of the catchment area 
(R?>0.89). This relationship was then adopted to simplify the classical Snyder method to determine UHs. 
Moreover, several parameters of the Snyder method were calibrated to the arid environment by matching 
the peak-flow, lag-time and three time-widths (75%, 50% and 30% of the peak-flow) of 10-min and 1-h 
UHs with physical characteristics of the watersheds. All developed relationships were validated with 
independent rainfall and wadi-flow events. Results indicate that the calibrated parameters in these arid 
watersheds are quite distinct from those suggested for other regions of the world. A marked difference was 
found between the 10-min UHs estimated by the S-hydrograph method and the deconvolution method. 
Therefore, it is concluded that a method depends on natural hydro-meteorological conditions would be 
more practical in arid region. The proposed methodology can be used for water resources management in 
arid regions having similar climate and geographical settings. 
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1 Introduction 


Floods are perhaps the most concern natural phenomenon in arid regions, which proved to have 
both positive and negative impacts. Moderate flood, which spread over a large area and been 
regulated by wide-opened channels, successive dikes and retention basins, can significantly 
replenish the groundwater aquifers (Abdalla et al., 2013, 2014; Mayaud et al., 2019). Young et al. 
(2004) discovered freshwater lenses in the Central Oman shallow aquifers claimed to be formed by 
recharge water from infrequent cyclonic rainfall. In contrast, flash floods are rare but have been a 
potentially catastrophic hazard in many parts of the world (Badrzadeh et al., 2015; Ahmadalipour 
and Moradkhani, 2019) including the arid regions. For example, Cyclone Gonu in 2007 is 
considered to be the worst natural disaster in recorded history of Oman, which caused 
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approximately 4.6 billion USD damage and 69 deaths. 

Flash floods are affected by various hydro-metrological characteristics, most notably triggered 
by extremely heavy rainfall, but also intensify by steep slopes, poor vegetation cover, soils with 
low permeability and ephemeral streams (locally known as wadi) attributed to long dry periods 
(Ghoneim and Foody, 2013; Parsons, 2018). Damages caused by the flash floods exaggerate as 
population growth and asset progressively expands in flood-prone areas (Hering et al., 2004; 
Kazama et al., 2009; Merz et al., 2010). Most of the areas categorized into arid climate zones in 
Middle East region receive annual total rainfall less than 130 mm with occasional high intensity 
and short-duration storms, which can generate river flows rise from zero to peak within a few 
minutes up to a few hours (Greenbaum et al., 1998). 

In addition to the difficulty with lack of gauging stations, river-flow simulation in arid regions 
present unique challenge due to rare occurrence of rainfall events and their greater spatial 
variability (Pilgrim et al., 1988). In response to these issues, hydro-metrological data network has 
been expanded and updated in various countries (Sen, 2008). As a result, many of the large 
watersheds in Oman have been equipped with river-gauge stations at their outlets, which provides 
an important opportunity to understand complex rainfall-runoff relationship in arid regions and will 
eventually enhance the clarity of water management and adaptation measures. Owing to population 
growth and industrialization, urbanization gradually expands from coastal valley areas to the 
upstream catchment areas in Oman. Planning and designing of various infrastructures such as roads, 
bridges, culvert and recharge dams warrant wadi-flow data, which are not yet available for many 
small and medium scale watersheds. 

Traditional runoff simulation methods such as Sherman's unit hydrograph (UH) technique 
(Sherman, 1932) and Tank model (Sugawara, 1995) use concurrent rainfall and wadi-flow data to 
develop complete flood hydrographs at the catchment outlet resulted from arbitrary storms. In the 
absence of river-flow and rainfall data in ungauged catchments, synthetic UH method is widely 
used for hydrologic/hydraulic designs (Bhuyan et al., 2015; Tomirotti and Mignosa, 2017). As 
explained by Chow et al. (1988), there are three types of synthetic unit hydrographs (UHs): (1) 
those used geomorphological properties of the watershed to develop theoretical or empirical 
relationships for the UH characteristics such as the peak-flow rate and timing etc. (Snyder, 1938; 
Gray, 1961); (2) those based on dimensionless UH (Soil Conservation Service, 1972); and (3) those 
based on models of watershed storage (Clark, 1943). 

Synthetic UH method does not depend on hydrological observations and is widely used in 
ungauged watersheds. However, parameters appeared in empirical relations are subjected to a great 
level of uncertainty when the original equations are applied with radically different hydro- 
meteorological conditions. For example, Snyder (1938) developed a set of empirical equations for 
the time-lag and peak-flow of his proposed synthetic UH by analyzing runoff characteristics in a 
large number of drainage-areas (26-26,000 km?) in Appalachian Highland in eastern United State. 
Since then, many studies have been conducted to regionalize model parameters in many parts of 
the world. Borga et al. (1991) used the Snyder method for developing a real-time flood forecasting 
system in Italy. Usul and Tezcan (1995) developed UHs from runoff data for four catchment areas 
in Turkey. The characteristics of these UHs were then matched with the Snyder's synthetic UH to 
obtain a set of parameters best-fit with local site conditions. Jena and Tiwari (2006) applied a 
similar methodology in two medium sized agricultural watersheds and a set of model parameters 
applicable for West Bengal state in India were proposed. El Hassan et al. (2013) used the Snyder 
method embedded in the HEC-HMS software to simulate rainfall-runoff relationship of a semi- 
urbanized catchment in Texas of the United States. 

The main objective of this study was to develop a methodology that is easily adaptable for 
designing hydraulic structures in ungauged catchment areas. A set of type curves were presented to 
develop the UH based on three simple but important catchment characteristics: watershed area, 
length of the main wadi and length to the centroid of the watershed along the main wadi from its 
outlet. This relationship was developed based on Snyder's UH method but was easily 
operationalizable for practical designs in ungauged catchment areas. Most of the past studies 
selected coarse time durations (e.g., 1-h or 2-h) when developing UHs due to the absence of rainfall 
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and wadi-flow in finer time resolutions. UHs with coarse time resolutions might not be accurately 
convertible to finer time durations to match with the intense and short-duration rainfall events 
common in arid regions. This study used wadi-flow data accumulated at 10-min intervals at the 
outlets of 12 watersheds with catchment area ranges from 45 to 1660 km? in Oman. Concurrent 
rainfall data at 1-min interval collected using tipping bucket rain gauges at 37 stations were also 
obtained. These data at finer time resolutions provided a unique opportunity to understand complex 
hydrological response in arid regions (1-h and 10-min UHs). 


2 Study area 


Oman is located at the southeast strip of the Arabian Peninsula (Fig. 1). Densely populated urban 
cities, such as Sohar, Muscat (the capital city), Sur and Salalah are located along the 3165 km long 
coastal plain. A series of rugged mountains (with the highest peak of 3075 m above sea level) and 
the Ophiolite hills run parallel to most of the coastline. The interior region accounts for nearly 82% 
of the country and mainly consists of sandy deserts. The climate in Oman varies mainly from semi- 
arid to hyper-arid. According to the records of 49 gauge stations over the 1985—2004 period, total 
annual rainfall ranges from 31 mm in the interior desert to 330 mm in the northern Oman mountain 
region, with a station average of 129 mm/a. The number of wet days (rainfall more than 1 mm) 
over the same period varies from 2 to 32 d per year, with a station average of 12 d (Gunawardhana 
et al., 2018). 
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Fig. 1 Delineated catchments and locations of rain gauge stations (a) and Thiessen Polygons for Aday catchment (Б) 


AlSarmi and Washington (2013) analyzed rainfall records at 23 stations covering 6 countries in 
the Arabian Peninsula including Oman for the period 1970—2008 and showed that the number of 
wet days with rainfall more than 10 mm significantly decreased. Gunawardhana and Al-Rawas 
(2016) studied rainfall records during the period 1990—2010 in three major coastal cities in Oman 
and found that the contribution from extreme rainfall to the annual total steadily increased. 
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Moreover, Gunawardhana et al. (2018) predicted enhanced tropical cyclone activities in the 
northeast coastal region during 2040—2059 and 2080—2099 periods. This study also predicts 
decreases in number of rainy days and increases in rainfall intensity in many regions of the country. 
These research findings highlight the necessity of rainfall-runoff modeling in this region to cope 
with potential changes in rainfall regime in the future. 


3 Methodology 


The research methodology consists of three major steps: (1) delineation of the catchment areas and 
river-network using digital elevation model (DEM) data and extraction of watershed 
characteristics; (2) derivation of UHs at different time durations from observed river-flow and 
rainfall records; and (3) matching the physical characteristics of the watersheds with features of the 
corresponding UHs to develop hydro-morphological relations. 

In the first stage, paired rainfall and river-flow data were collected to cover 12 major watersheds 
(Fig. 1). Observed river-flow at 5-min interval and tipping bucket rainfall measurements at 1-min 
interval were obtained from the Ministry of Regional Municipalities and Water Resources of Oman. 
Selected hydrograph cover a large extent of river-flow rates including the highest ever recorded 
during hurricane Gonu in 2007. Key catchment characteristics such as the watershed area, length of 
the main river and length to the centroid of the watershed along the main river from its outlet were 
derived from 40-m resolution DEM data and HEC-GeoHMS (Geospatial Hydrologic Modeling 
Extension) tool. The selected watersheds (Table 1) have a wide range of areas (45-1660 km’) and 
river lengths (13—92 km), which is a key need for generalizing results from small to large ungauged 
basins. In general, shape of the catchment area affects the magnitude and timing of the peak-flow 
(Nouh, 1990). The whole study area, which includes 12 major watersheds were then divided into 515 
sub-catchments with various shapes and sizes. Basin shape factor, which was computed by dividing 
square of the river length by watershed area (Al-Rawas and Valeo, 2010), ranges from 0.5 (very 
circular) to 19.5 (extremely elongated). Similarly, sub-catchment areas extends over a wide range of 
0.01 to 42.00 km”. 


Table 1 List of stations and delineated catchment characteristics 


Watershed Wadi gauge Area (Кт?) Wadi length (km) Length to the centroid (km) 
Aday FB505467AD 794.2 54.6 25.2 
Afi EL895407AD 313.0 34.2 16.3 
Al Khawd FB104840AD 1660.2 91.9 39.6 
Al Uqq FA160968AD 132.1 22.5 8.2 
Arabiyin GA057335AD 303.8 37.7 18.9 
Hayfadh FA777631AD 122.2 21.7 7.4 
Lansab FA395799AD 44.7 13.3 8.5 
Manzariya FA585595AD 206.4 28.3 9.8 
Mayh FA596055AD 222.0 34.0 15.5 
Miglas FA877343AD 556.8 55:5 25.7 
Awabi EL574613AD 254.4 35.4 14.4 
Jannah FA580672AD 130.8 29.9 18.9 


In the second stage, river-flow observations were aggregated to 10-min and 1-h durations when 
generating river-flow hydrographs. Similarly, 10-min and 1-h total rainfalls were calculated from 
tipping bucket rain gauge records to match with river-flow durations. When there are more than 
one rain gauge station correspond to a catchment area, the Thiessen Polygon method was used for 
averaging the rainfall time series. For example, the Thiessen Polygons developed for the Aday 
catchment area are shown in Figure 1b. There are no perennial river-flow in these watersheds. In 
fact, most of the watersheds in arid regions have ephemeral-type rivers that produce water flow 
temporally during the rainfall and dried up immediately due to extremely high evaporation demand 
and infiltration rate. Therefore, base-flow is essentially negligible and excess rainfall volume is 
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equal to the integrated area under the hydrograph. The reverse process of the discrete UH 
convolution technique, known as the deconvolution method (Chow et al., 1988), was applied to 
derive a UH from river-flow and excess rainfall. Excess rainfall was calculated by the constant 
fraction loss model (McCune, 1998; Dingman, 2015), where a certain percentage of rainfall become 
runoff. Corresponding mathematical equation can be expressed as: 


i(t)=C, p(t) 
[о()а, (1) 
C, == 
A f p(t)dt 
where i(t) is the excess rainfall rate (m/h); p(t) is the observed rainfall rate (m/h); C, is the fraction 
of rainfall converts to runoff; Арда is the rainfall volume (m°); A is the watershed area (m7); and 
| O(t)dt is the cumulative wadi-flow volume (m°). 
In the third stage of this research, the basics of the Snyder's empirical relationship was followed 
to relate UH features with the physical characteristics of the watershed. According to Snyder 


(1938), lag-time (tp), which is defined as the time from the center of mass of excess rainfall to the 
peak-flow rate, is given by: 


1 =C,(LxL,)’, (2) 


where L is the main wadi length from the most upstream point of the watershed to the catchment 
outlet (km); Le is the distance from the catchment outlet to the point along the main wadi which is 
the closest to the centroid of the watershed (km); and C; and n are parameters which depends on 
particular catchment characteristics. Snyder (1938) calibrated n=0.3 by analyzing a large number 
of hydrographs in Appalachian Highland of eastern United State. Since then, many studies used the 
same value without proper calibration into local conditions including in arid regions (Sen, 2008; 
Sudhakar et al., 2015). In this study, parameters п and С, were calibrated by matching the properties 
of the derived UHs. Also, according to Snyder (1938), the peak-flow rate is given as follows: 
(2.78C,, x A) 
ЕИ 

where Q, is the peak-flow rate (т2/); Су is a parameter which depends on the storage and retention 
characteristics of the watershed; and A is the catchment area (km7). In this study, a relationship 
between L and L. was developed by considering the 515 sub-catchments delineated in the first step. 
Because С, апа С, depend on particular catchment characteristics, they were combined together to 
form a single parameter (C in Eq. 4). It was later calibrated by matching the properties of the 
derived UHs. Accordingly, the relationship between L and Le will simplify the Equation 3, where 
О» can be estimated in terms of two physical features of the watershed: L and A as in Equation 4. 


_CxA 


Q, (L)" Ч 


where the exponent m can be calculated based on the relationship between L апа Le. 


(3) 


(4) 


4 Results and discussion 


4.1 Extraction of catchment characteristics from DEM data 


Until recently, basic catchment characteristics, such as the watershed area and wadi-length, were 
derived from topographic map sheets manually. Some watershed parameters, such as the slope and 
the centroid of the area, are not readily available through manual delineation, but must be 
calculated. It therefore is difficult to estimate them with sufficient accuracy. As DEM data and 
automatic catchment delineation become widely popular, reliability of the parameter estimations 
significantly improved. However, it requires certain experience and sophisticated software, which 
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may not be available yet among designers and practitioners at local government levels. Therefore, 
this study was planned to reduce the number of parameters required to estimate flood hydrograph, 
and presents it based on certain parameters, which are easily achievable. Accordingly, a relationship 
between L and L. was developed. Out of the 515 sub-catchments delineated, half of them were used 
to develop the relationship (Fig. 2a) and rest were used for model validation (Fig. 2b). This 
relationship was further validated with the information delineated from 12 major watersheds (Fig. 
2c). All three relationships included various catchment areas with different shapes, and produced 
higher coefficients of determination (R?>0.89). This relationship therefore can be successfully 
applied in Equation 3 to simplify the estimations of Q, (Eq. 4). 
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Fig. 2 Calibrated and validated relationship between L (the main wadi length from the most upstream point of the 
watershed to the catchment outlet) and Le (the distance from the catchment outlet to the point along the main wadi 
which is the closest to the centroid of the watershed) 


4.2 Derivation of UHs at different time durations 


Eight of the 12 major watersheds studied were used for developing the relationships to calculate 
lag-time and peak-flow rate of the UH (Eqs. 2 and 4, respectively) and remaining four watersheds 
were used for verifying those relationships. The observed wadi-flow rates at 10-min and 1-h 
intervals along with the calculated excess rainfall (Eq. 1) at same durations were used with the 
deconvolution method to estimate the corresponding UHs. In verification stage, above UHs were 
applied with different rainfall events to calculate resultant wadi-flow. Figure 3 presents the 
observed and simulated wadi-flows for calibration and verification of 1-h UHs. The UHs were 
calibrated using short duration rainfall events, yielding a well-defined single peaked hydrograph of 
short time base, and verified using the events with the highest observed peak-flow. Four of the eight 
major watersheds were verified with observed wadi-flows occurred during Cyclone Gonu in 2007. 
When there are more than two hydrographs available for a catchment area, they were all used for 
verifying the calibrated UHs. The coefficients of determination (R?>0.89) evident for a good 
agreement between observed and simulated wadi-flows. Some level of differences between 
observed and simulated wadi-flows can be attributed not only to the limitations of the model itself 
but also to the various issues such as observational errors of wadi-flows and rainfalls, and 
inadequacy of the meteorological stations to capture the spatial distribution of the rainfall events. 
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Fig. З Observed and simulated wadi-flow using 1-h UH (unit hydrograph) 


When a UH of a given excess-rainfall duration is available, the principle of superposition (S- 
hydrograph method) can be used to derive the UHs of any required duration (Chow et al., 1988). 
Accordingly, 1-h UHs developed in the previous stage were used to calculate 10-min UHs. Also, 
another set of 10-min UHs were derived using 10-min rainfall and wadi-flow data applied with the 
deconvolution method. Figure 4 depicts a marked difference of the 10-min UHs developed by these 
two methods for Aday watershed. Similar phenomenon in different magnitudes were observed in 
all watersheds where the 10-min UH developed by the 10-min rainfall and wadi-flow data produced 
a greater peak-flow. The reason for this difference can be attributed to the assumptions made in 
deriving UHs. For example, it is assumed that the excess rainfall is uniformly distributed over the 
entire watershed at a constant rate for the selected duration. This assumption, under natural 
conditions in arid regions, can rarely be satisfied. Arid regions experience occasional rainfall events 
with intensity varies significantly over time (Greenbaum et al., 1998). Therefore, such difference 
between theoretical S-hydrograph method and the deconvolution method with real data can be 
expected. Figure 5 shows the peak-flow simulated by these two methods and their match with the 
observed 10-min averaged peak-flow. It is noted that the UH derived from the S-hydrograph 
method underestimated the peak-flow up to 30% while the deconvolution method with short-term 
real data produced promising results. It is also clear that this difference between the observed and 
simulated peak-flows by the S-hydrograph method increases as the catchment area increase. These 
results therefore demand appropriate modifications for existing methods used for wadi-flow 
estimation in arid regions. 
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Fig. 4 Comparison of 10-min UHs developed by two methods for Aday catchment 
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Fig.5 Observed and simulated 10-min averaged peak-flows by two methods for all watersheds 


4.3 Relationship of UH features with physical characteristics of watersheds 


4.3.1 Peak-flow 


According to Snyder (1938), peak-flow (Qp) can be expressed as a function of three physical 
characteristics of a watershed: area (A), length of the main wadi (L) and length to the centroid of 
the watershed (Lc) (Eq. 3). Based on the relationship developed between L and Le in Figure 2a, 
Equation 4 can be used to estimate Q, as a function of A and L. Accordingly, the peak-flows taken 
from the derived UHs and delineated catchment characteristics (Table 1) were used with the least- 
square method to obtained parameters C and m. Figure 6a depicts a comparison of the peak-flows 
taken from the derived UHs and the ones derived with estimated C and m. Because the objective of 
this study is to develop a methodology that is easily adaptable in designs, it is important to keep 
calibrated parameters unified for different time durations. Accordingly, the calibrated m value of 
0.7 was applied with the peak-flows from the derived 10-min UHs. Figure 6b shows that the peak- 
flow matches with the A/L°7 ratio reasonably well (R?=0.98). Hence, т=0.7 can be used in 
estimating peak-flow of UHs at 1-h апа 10-min durations, and it might be applicable for other 
durations as well. The differences between the 1-һ and 10-min UHs for various features, such as 
the peak-flow and lag-time, can be accounted by markedly different gradients of the two graphs 
(1.89 and 3.01 m°/(s-km!) for 1-h and 10-min UHs, respectively). 
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Fig. 6 Relationship between peak-flow (Оз)апа physical characteristics of the catchment derived using the 
optimal m value (0.7) for 1-h and 10-min UHs (unit hydrographs) 


4.3.2 Lag-time 


The lag-time (tp) is defined as the time from the center of mass of excess rainfall to the peak-flow 
rate, which according to Equation 2 can be expressed as a function of the length of the main wadi 
(L) and length to the centroid of the watershed from its outlet (Le). Theoretically, the lag-time of 
the 10-min UH should be smaller than that of the 1-h UH (Fig. 4). When the п parameter of the 
Equation 2 is calibrated (n=0.35) and fixed for any time duration, remaining parameter С, should 
be varied to account the difference between 10-min and 1-h UHs. Figure 7 shows the calibrated С, 
for estimating lag-time for the 1-h and 10-тіп UHs. The relationship between t, and (LxL.)°*> for 
the 1-h UH produce only a moderate agreement with a coefficient of determination equals to 0.74. 
This is because the lag-time for 1-h UHs can only be estimated with a +0.5 h accuracy. When actual 
lag-time of different catchments varies within a smaller range (1.0-3.5 h as in Fig. 7), such 
uncertainty in estimating tp can be significant. Due to the finer time resolution, the relationship 
between tp and (LxL,)°* for the 10-min UH produces good agreement with a coefficient of 
determination equals to 0.98. 
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Ее. 7 Calibration of С, (Equation 2) for estimating lag-time for the 1-h and 10-min UHs 


4.3.3 Time-widths 


Sections 4.3.1 and 4.3.2 present quantitative analyses to estimate the peak-flow and lag-time. 
However, in order to determine the shape of the UH, few more ordinates on the UH must be known. 
Snyder (1938) expressed the time-width of the UH at a discharge equal to a certain percent of the 
peak-flow as a function of peak-flow rate per unit watershed area. It was further assumed that one- 
third of this width is distributed before the UH peak time and two-third after the peak (Chow et al., 
1988; Angelidis et al., 2010). Accordingly, three relationships were developed for the time-width 
at 75%, 50% and 30% of the peak-flow based on the derived UHs (Fig. 8). It is clear that the time- 
width is inversely proportional to the peak-flow per unit watershed area. All six graphs produced 
relatively acceptable relationships with determination coefficient more than 0.87. These empirical 
equations provide six more points on the UH, which brings altogether eight ordinates to determine 


202101.00055v1 


chinaXiv 


ChinaXivA ERAT 


JOURNAL OF ARID LAND 


the probable shape of the UH. 
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Fig. 8 Time width at 75% (№5), 50% (Жо) and 30% (W30) of the О for 1-h and 10-min UHs. A, catchment area. 


Rainfall-runoff hydrographs in arid regions can be quite distinct from those in humid regions. 
Attributed to the steep slope, flow is faster in arid regions resulting in shorter rising limb and greater 
peak-flow. The total discharge significantly decreases within short periods of time due to high rates 
of evaporation and transmission loss. Wadi-beds in these regions are generally composed of alluvial 
sediments dominated by highly permeable coarser grain sizes. Therefore, 1:2 width-ratio as 
suggested by Snyder (1938) may not be applicable in arid regions. A series of UHs was developed 
for 1-h and 10-min durations for all study areas based on empirical relationships developed in 
previous steps (Figs. 2, 6, 7 and 8). When compare these UHs with ones derived from observed 
wadi-flow, it was determined that the 3:4 width-ratio produced better results than the 1:2 width- 
ratio. For example, Figure 9 shows a comparison of 10-min UHs developed for Afi watershed with 
two different width-ratios. It is clear that the 3:4 width-ratio formed shorter rising limb than the 1:2 
width-ratio and matches better with the derived UH from observed wadi-flow. 
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Fig. 9 Comparison between developed 10-min UHs for Afi watershed with two different width-ratios 


The applicability of the proposed empirical equations to develop the UHs at 1-h and 10-min 
durations was tested by simulating wadi-flows in four more watersheds. These four watersheds are 
different to those used for calibrating and verifying UHs in Figure 3, whose catchment 
characteristics are shown in Table 1. Figure 10 shows the observed and simulated hydrographs 
based on the empirical equations proposed. It is clear that the simulations by the 1-h UHs agree 
better with the observations than the results from the 10-min UHs. When simulating wadi-flows at 
finer time resolutions, uncertainty arises due to temporal and spatial variations of rainfall events, 
especially in arid regions. However, time averaged 1-h wadi-flows smooth the abrupt variations at 
10-min interval and therefore produced better agreements. All simulations reasonably capture 
different features of the hydrographs, such as the general trend, peak-flow, time-lag etc., when 
match with the observations. It therefore indicates that the developed empirical relationships can 
be applied in deriving UHs particularly for arid regions in Oman. These results can further be 
improved by applying them in different countries with similar climate and geographical settings. 


5 Conclusions 


Rainfall-runoff relationship in arid regions is unique and proper understanding is challenging. This 
study aimed in developing a simple but reliable procedure for estimating UH at reasonably finer 
time resolutions (10-min and 1-h) which can be easily adaptable by practitioners at sub-catchment 
levels, especially when it is ungauged. For example, wadi length can easily be estimated by 
automatic delineation from DEM data or even by manual methods. However, calculation of the 
catchment centroid is rather difficult, especially when paper based topographic maps are used. 
Results of this study indicate that the strong relationship found between the wadi length and the 
length to the centroid of the catchment area can be adopted to simplify the existing methods to 
determine UHs. Moreover, it was found that the principle of superposition with the assumption of 
uniform excess rainfall distribution would not be applicable in arid regions. A marked difference 
was found between the 10-min UHs estimated by the S-hydrograph method and deconvolution 
method. Therefore, it is concluded that a method depends on natural hydro-meteorological 
conditions would be more practical in arid region. 

A set of empirical equations were developed for estimating peak-flow, lag-time and time-widths 
at 75%, 50% and 30% of the peak-flow for 10-min and 1-h UHs by analyzing wadi-flows in eight 
watersheds. Developed relationships were successfully validated with observed wadi-flows in four 
independent watersheds. When calibrated parameters are concerned, developed empirical 
relationships are quite distinct from those suggested for other regions. For example, the 3:4 width- 
ratio was found to be more suitable in arid regions than 1:2 width-ratio suggested in numerous past 
studies. The proposed methodology is expected to be useful for water resources management in arid 
regions, especially for planning countermeasures to cope with future climate change. 
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Mayh: Verification of 1-h UH (2 Nov 2011) Mayh: Verification of 10-min UH (2 Nov 2011) 
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Fig. 10 Comparison of observed and simulated wadi-flows from proposed empirical equations 
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